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Abstract. Recent work has hypothesized that tropical cy- 
clones in the deep Atlantic and eastern Pacific basins de- 
velop from within the cyclonic Kelvin cat’s eye of a trop- 
ical easterly wave critical layer located equatorward of the 
easterly jet axis. The cyclonic critical layer is thought to be 
important to tropical cyclogenesis because its cat’s eye pro- 
vides (i) a region of cyclonic vorticity and weak deforma- 
tion by the resolved flow, (ii) containment of moisture en- 
trained by the developing flow and/or lofted by deep convec- 
tion therein, (iii) confinement of mesoscale vortex aggrega- 
tion, (iv) a predominantly convective type of heating profile, 
and (v) maintenance or enhancement of the parent wave until 
the developing proto-vortex becomes a self-sustaining entity 
and emerges from the wave as a tropical depression. This 
genesis sequence and the overarching framework for describ- 
ing how such hybrid wave-vortex structures become tropical 
depressions/storms is likened to the development of a marsu- 
pial infant in its mother’s pouch, and for this reason has been 
dubbed the “marsupial paradigm”. 

Here we conduct the first multi-scale test of the marsupial 
paradigm in an idealized setting by revisiting the Kurihara 
and Tuleya problem examining the transformation of an east- 
erly wave-like disturbance into a tropical storm vortex using 
the WRF model. An analysis of the evolving winds, equiv- 
alent potential temperature, and relative vertical vorticity is 
presented from coarse (28 km), intermediate (9 km) and high 
resolution (3.1 km) simulations. The results are found to sup- 
port key elements of the marsupial paradigm by demonstrat- 
ing the existence of a rotationally dominant region with min- 
imal strain/shear deformation near the center of the critical 
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layer pouch that contains strong cyclonic vorticity and high 
saturation fraction. This localized region within the pouch 
serves as the “attractor” for an upscale “bottom up” develop- 
ment process while the wave pouch and proto-vortex move 
together. 

Implications of these findings are discussed in relation to 
an upcoming field experiment for the most active period of 
the Atlantic hurricane season in 2010 that is to be conducted 
collaboratively between the National Oceanic and Atmo- 
spheric Administration (NOAA), the National Science Foun- 
dation (NSF), and the National Aeronautics and Space Ad- 
minstration (NASA). 


1 Introduction 

The genesis of tropical cyclones has been regarded, until very 
recently, as one of the “great mysteries of the tropical atmo- 
sphere” (Emanuel, 2005). Like a steep unconquered peak in 
the Himalaya, for which a tractable route of ascent remains 
to be discovered, the genesis problem presents a number of 
scientific and technical challenges from observational, theo- 
retical and numerical directions. The multi-scale nature of 
the problem has been arguably the most daunting crux of the 
climb. We believe that a successful summit bid requires a 
quantum leap in dynamical understanding of the Lagrangian 
nature of the genesis process and that such insight will open 
new routes to the top. 

Building on previous work examining tropical easterly 
wave lifecycles (e.g. Thorncroft and Hoskins, 1994) and up- 
scale organizational processes by vortical hot towers (VHTs) 
within such favorable environments (Montgomery et al., 
2006), recent work by Dunkerton et al. (2009; hereafter 
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DM W09) proposed a new paradigm for tropical cyclogenesis 
that recognizes the intrinsic multiscale nature of the problem 
from synoptic, sub-synoptic, mesoscale and cloud scales. 
Using three independent datasets, ECMWF Reanalysis data, 
TRMM 3B42 3-h precipitation and best track data from 
the National Hurricane Center (NHC), the Kelvin cat’s eye 
within the critical layer of a tropical easterly wave was hy- 
pothesized to be important to tropical storm formation be- 
cause: HI) wave breaking or roll-up of cyclonic vorticity 
and lower-tropospheric moisture near the critical surface in 
the lower troposphere provides a favorable environment for 
the aggregation of vorticity seedlings for TC formation; (H2) 
the cat’s eye is a region of approximately closed circulation, 
where air is repeatedly moistened by deep moist convection 
and protected to some degree from dry air intrusion; (H3) 
the parent wave is maintained and possibly enhanced by dia- 
batically amplified mesoscale vortices within the wave. This 
genesis sequence and the overarching framework for describ- 
ing how such hybrid wave-vortex structures become tropi- 
cal depressions is likened to the development of a marsu- 
pial infant in its mother’s pouch wherein the juvenile proto- 
vortex is carried along by the mother parent wave until it 
is strengthened into a self-sustaining entity. A survey of 
55 named storms in the Atlantic and eastern Pacific sectors 
during August-September 1998-2001 was shown to support 
this so-called “marsupial paradigm’'. Tropical cyclogenesis 
tended to occur near the intersection of the trough axis and 
the critical surface of the wave, the nominal center of the 
cat’s eye. 

The marsupial paradigm provides a useful road map for 
exploration of synoptic-mesoscale linkages essential to trop- 
ical cyclogenesis. In this paper we begin a systematic test 
of the marsupial paradigm by revisiting the canonical prob- 
lem of the transformation of an idealized easterly wave-like 
disturbance into a tropical storm vortex. Both coarse res- 
olution (cloud parameterization) and high resolution (cloud 
representing) numerical experiments are conducted, and the 
high-resolution experiment includes a region of intermedi- 
ate horizontal resolution in order to compare, side-by-side, 
two otherwise identical storms with near cloud-resolving nu- 
merics. These experiments are a first step towards closing 
the gap between observational I ly resolvable scales (synop- 
tic, meso-alpha) and infrequently sampled cloud-system and 
cloud-scale process (meso-beta and meso-gamma, including 
VHTs) 


□ Model and eEJjeriment design 

□ .1 WDD model 

The numerical model used is the Weather Research and Fore- 
casting (WRF) model, developed and maintained by NCAR 
(Skamarock et al., 2005). The WRF model is a fully com- 
pressible, Euler nonhydrostatic model with a run-time hy- 
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drostatic option available. It adopts Arakawa C-grid stag- 
gering to facilitate the natural discretization of the continu- 
ity and tracer equations. Time-split integration uses a 3rd- 
order Runge-Kutta scheme with smaller steps for acoustic 
and gravity wave modes, and 2nd to 6th order advection op- 
tions are used in the horizontal and vertical discretizations 
in the form of flux divergence. The model uses a terrain- 
following, hydrostatic-pressure vertical coordinate, which is 
also called a “mass vertical coordinate”. 

Our WRF implementation employs the nonhydrostatic op- 
tion and possesses 33 vertical levels with 23 of them below 
500 hPa. Symmetric lateral boundary conditions are imposed 
at the north-south boundaries, and periodic boundary condi- 
tions are used in the west-east direction. Standard “physics” 
packages are employed for the parameterization of deep 
and shallow convection (Betts-Miller- anjic scheme), micro- 
physics (Kessler scheme), PBL ( onsei University scheme), 
longwave radiation (RRTM scheme) and shortwave radiation 
(Dudhia scheme). 

□ □□ primental design 

□.□•l Dasic state (Ibw 

Our objective is to test the hypothesized role of the critical 
surface and its associated finite-amplitude cat’s eye flow in 
the ensuing multi-scale convective evolution within and near 
the easterly wave’s pouch. To help provide a benchmark 
for future studies, the now classical tropical cyclone gene- 
sis experiment performed by Kurihara and Tuleya (1981 ; re- 
ferred to hereafter as KT81) will be re-conducted here and 
interpreted with new eyes using the marsupial paradigm of 
DMW09. 

We adopt the same basic state flow as in KT81 . The initial 
basic flow is defined as follows: 

u(<p, cr)=Micos(or^)tanh(^ a ) + Mo^^ 

U = 0 

where <7 = ^-, P* is the surface pressure, u and v denote the 
zonal and meridional flow of the base state, respectively, 0 
denotes latitude, 0o 1 8° N is the center latitude of the zonal 
channel, and D 2.5°, u\ 3.1 5 ms -1 and uo —5 ms -1 . 
This idealized zonal flow is consistent with the observed 
zonal flow during Phase III of GATE near 80° W (and con- 
sistent also with recent NCAR re-analyses for this region be- 
tween 1968-1998 as reported in Fig. I of Hall et al., 2006). 
The main difference between this model configuration and 
that of KT81 is that the channel width is extended to 50° lat- 
itude in order to reduce the influence of the lateral boundary. 
Following KT81, the mean surface pressure, temperature and 
geopotential are obtained from the basic flow by inverting the 
geostrophic and hydrostatic balance equations in Cartesian 
coordinates (see Appendix A for details). 

As shown in Fig. 1, an easterly mean flow is prevalent 
throughout the model domain, and a broad easterly jet is 
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□ig. 1. Latitude-height distributions of the idealized basic zonal 
flow (top, ms -1 ), the corresponding relative vorticity (middle, 
10“ 6 s -1 ) and the meridional gradient of the relative vorticity (bot- 
tom, 10“ 11 m _1 s _1 ). The ordinate is or and the abscissa is latitude. 

present in the northern part of the domain. The horizontal 
shear of the zonal flow is largely confined within ^400 km 
on both sides of the center latitude (18° N). North of the cen- 
ter latitude, the basic flow has a westerly vertical shear; south 
of the center latitude, the basic flow has an easterly vertical 
shear of similar magnitude. Along the center latitude, the 
basic zonal flow is -5 ms -1 and is constant with height (no 
vertical shear). Since the meridional gradient of the relative 
vorticity changes sign at the center latitude in both the upper 
and the lower troposphere, it satisfies the (Rayleigh) neces- 
sary condition for barotropic instability (Holton, 1992). For 
simplicity, the calculations here are carried out on an f-plane 
(corresponding to 18° N). (When the beta effect is included, 
the location where the meridional gradient of the absolute 
vorticity changes sign is slightly north of the center latitude). 

□ (Unstable modes and initial disturbances 

To determine the most unstable (dry) inviscid mode for this 
zonal mean flow, small amplitude random noise to the ini- 
tial wind field was introduced, and the model was integrated 
forward in time without moist processes or surface friction. 
As is usual, with a periodic boundary condition in the east- 
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west direction, the zonal scale of the domain constrains the 
allowed wavelength of the unstable modes (the zonal scale 
of the domain is an integer multiple of the wavelength of any 
wave it resolves). So as to permit the excitation of a suffi- 
ciently broad spectrum the zonal domain is set to 100° longi- 
tude. To reduce computation time, the horizontal grid spac- 
ing was set to 56 km. The initial disturbance was produced 
by adding random wind noise at each grid point with maxi- 
mum magnitude about 0.25 m s _1 (equivalent to random rel- 
ative vertical vorticity disturbances of order 10 -5 s _1 ). The 
magnitude of the wind noise, A, is taken to be independent 
of height, but decays exponentially away from the center lat- 
itude and is confined within 8.5° on both sides of the center 
latitude (Fig. 2a) according to the formula: 

_ { Ao[exp(-f(f-) 2 )-exp(-|)] if lf;l<I 
A — 1 , . y 

I 0 otherwise 

where L y i 952 km (8.5°), Ao 0.25 ms -1 , and y is zero at 
the center latitude (18° N), where the meridional gradient of 
the relative vorticity changes sign. 

After 45 days of integration two wave groups dominate the 
response. The first group is confined in the lower troposphere 
(below 500 mb) and the maximum perturbation winds reside 
near the lower surface (Fig. 2b). The second group of waves 
is confined in the upper troposphere (not shown), with a sim- 
ilar magnitude and a longer wavelength (~3 00km). By 35 
days the perturbation meridional velocity magnitude at the 
surface increases to approximately 1 ms -1 , and the wave- 
length is approximately 25° longitude (about 2800 km). As 
shown in Fig. 3, the e-folding time of this lower-tropospheric 
mode is approximately 9.8 days, but it takes approximately 
30 days, including the incubation period, for the random per- 
turbations to develop into a coherent wave train pattern. As 
is evident from Fig. 2b, these waves tilt SW-NE against the 
prevailing latitudinal shear during their growth phase, sug- 
gesting a barotropic instability. After 45 days of integra- 
tion Figs. 2c and 3), the growth of the dominant mode in 
the lower troposphere has ceased and the mode has reached 
a saturation amplitude for the meridional wind of approxi- 
mately 1.8ms -1 . The numerical model conserves the cor- 
responding local extremum of vorticity reasonably well at 
subsequent times (not shown). From Fig. 2c we see that the 
flow evolution has transformed the initial vortex strip into a 
train of (cyclonic) Kelvin cat’s eyes with maximum cyclonic 
vorticity inside each closed “gyre-pouch”. 

Given the slow development of dry waves in the KT81 
basic state, we conclude that this basic flow supports only 
a marginally unstable mode. The zonal phase speed of 
this marginally unstable wave is found to be approximately 
—5 ms -1 or 4 degrees longitude per day westward. This is 
the speed of the zonal mean zonal flow near the center lati- 
tude, and thus the critical latitude of the most unstable mode 
(and emergent Kelvin cat’s eyes) is very close to the center 
latitude of the model domain ( 1 8° latitude). 
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□ig. □ Meridional wind disturbance (unit: m s -1 ) at sea surface from the idealized dry model simulation (see text for details) at Da Pay 1 
and [bfDay 35, and SQ-elative vorticity (unit: 10 -5 s -1 ) and streamlines at Day 45 in the co-moving frame of reference that moves at the 
zonal phase speed of the waves (— 5 m s _I ). 



□g. □ Time series of the surface maximum meridional wind (black; 
y-axis on the right) and its natural logarithm (blue; y-axis on the 
left) from the dry WRF simulation. The red straight line is a linear 
fit on a log-linear scale from Day 16 to Day 35, implying an e- 
folding time of approximately 9.8 days. 


As a check on whether the domain is large enough to select 
the most unstable mode for this zonal flow, we extended the 
domain zonally to 120° longitude, and the dominant wave 
response did not change. This suggests that the wave with 
2800 km wavelength (or 25° longitude) is indeed the most 
unstable wave for this basic flow, subject to the integer mul- 
tiple constraint: viz., wavenumber 4 or 5 in the 100 or 120 
deg domains, respectively. The upper tropospheric distur- 
bance is wavenumber 3 in the 1 00 deg domain and evidently 
decoupled from the lower tropospheric wave. 

It might be argued that having to wait 45 days for the 
emergence of a nascent gyre-pouch with maximum cyclonic 
relative vorticity of only 1.5xl0 -5 s -1 is far too long to be 
of practical relevance for tropical cyclogenesis in the West- 
ern Atlantic or other basins such as the Western Pacific for 
which the zonal mean flow might serve as a reasonable ap- 
proximation in certain flow regimes. However, these dry 
lifecycle experiments usefully illustrate how easily a ran- 
dom wind perturbation can project onto the intrinsic cat’s eye 
flow structures that are supported by these types of cyclonic 
shear flows. 
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For the remainder of this paper we adopt the same method- 
ology as KT81 and investigate the disturbance-to-storm 
transformation process by employing a small but finite am- 
plitude initial perturbation having approximately the same 
structural characteristics (horizontal wavelength and verti- 
cal decay scale) of the dominant unstable (dry) wave pattern 
found above, without the meridional phase tilt, and at ap- 
proximately three times the saturation amplitude for relative 
vorticity (^4.5 x 10 -5 s' 1 ). The corresponding moist exper- 
iments conducted herein therefore adopts the “dumpiness” 
of the cyclonic gyre-pouches from the dry experiment, but 
magnifies their value by approximately a factor of three. Al- 
though the genesis sequence presented below has been ver- 
ified to be essentially similar to its small-amplitude coun- 
terpart experiment, this scaling renders our set up virtually 
identical to that of KT81 and gives an arguably more realis- 
tic transformational time scale from a precursor disturbance. 

Following KT81, the zonal domain is set to 25° longitude 
and the initial wave disturbance is defined from the following 
expressions for the zonal and meridional wind perturbation in 
Cartesian coordinates: 


-9^A( < T)cos(^f)exp(-|(^) 2 ) 



0 otherwise 


( 3 ) 


^A(a) sin(2gf) exp(-f(^)^-exp(-f) 1 


= • if |£|<1 


( 4 ) 


0 otherwise 


where L y 952 km (8.5°), L x 2800 km (25°) and A{o) de- 
fines the vertical profile of the disturbance 


Following KT81, the surface pressure, three-dimensional 
temperature and geopotential fields are derived from the ini- 
tial wind field by solving the reverse balance equation in 
Cartesian-sigma coordinates (see Appendix A for details). 
Rather than adjust these initial fields to a steady-state fric- 
tional boundary layer (as in KT81) we let the WRF model 
perform this adjustment. The boundary layer adjustment oc- 
curs on a time scale on the order of l/v^/Cabs) (cf. Eliassen 
and Lystad, 19 1 ), where £ a bs * s absolute vertical vortic- 
ity of the zonal mean flow (f — dU /dv). For / corresponding 
to 1 8° N, the adjustment time is less than 39 h. 

Inclusion of moisture 

To include moisture, we specify a vertical profile of relative 
humidity (r), and then deduce the air temperature ( T ) and 
mixing ratio ( q ) from virtual temperature (7V) by iteration 
upon using the Magnus formulae: 

T= T i ( i+ e , v) 

lOOOfiiexpl^lT" — ko)/(T — £3)1 •* 1 >F o 
I000fi|exp|21.8 (T — 2 316)/(T- , 66 )| otherwise 

<<s= £4es /(p-e s ) < 6 ) 

q=(1sr 

where q s is the saturation mixing ratio, p is total pressure 
and e. v is the saturation vapor pressure in Pa. The constants 
are defined as follows: £o 2 3.15 K, e\ 0.6112, £2 l .61 J 
e 3 29.65 K, and £4 0.62198. 

Following KT81, we specify a horizontally homogeneous 
relative humidity field on constant pressure surfaces. This 
is a reasonable approximation to the observed moisture field 
over the West Atlantic. The vertical cross-section of mixing 
ratio is shown in Fig. 4 f. Because of the cyclonic (anticy- 
clonic) wave-like disturbance that is added to the zonal ba- 
sic state flow, a corresponding low (high) surface pressure is 
present near the domain center. From Eq. ( 6 b), this implies 
a slightly enhanced (reduced) mixing ratio near the center of 
the domain on constant pressure surfaces (Fig. 4f). 


1 , for cr > ao 

cos 2 (t^^), for cro>cr >gq — L„ (5) 
0 , for o < ao — L n 

where i>o 5ms _I , L„= 0.35, and oro=O.U Figure 4a-c 
shows plan views of the zonal, meridional and relative vor- 
ticity fields at the lowest sigma level over one wavelength. 
From Fig. 4d, e we see that the initial meridional and zonal 
wind disturbance is largely confined below the 500 hPa level 
and has maximum amplitude at the surface. As noted by 
KT8 1 (p 1 634), “The decrease of wave amplitude with height 
implies, under the assumption of a balance condition, lack of 
a cold core structure. In this respect the present disturbance 
resembles certain waves observed in the Pacific ... more 
closely than those in the Atlantic.” 


□ Wavetjtoueh tracking method 

In our previous diagnoses using coarse-resolution global 
model products (DMW09; Wang et al., 2009), Flovmoller 
diagrams of meridional wind perturbations or moisture field 
were used for wave tracking. So as to obtain consistent re- 
sults from the high-resolution simulations, the propagation 
speed of the wave/pouch is determined by tracking of the 
pouch centroid. The pouch centroid ( x , y) is defined as 

[*=* £*1 

pouch , s 

— 1 v- v j 

y = „T y ( - 

pouch 

where i denotes the grid points within the pouch and n is 
the total number of grid points within the pouch. We define 


A(or) = ^v 0 
2n 
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□ig. □ The left panels show the initial zonal wind (a: unit: ms -1 ), meridional wind (b. ms -1 ) and relative vorticity (c; unit: 10 -5 s -1 ) at 
the surface, and the right panels show the vertical cross section of the initial zonal wind perturbations S1Q meridional wind perturbations 3D 
and mixing ratio (f; unit: g/kg). Perturbations are defined as the difference between the total field and the basic state. 


the pouch by the closed stream function contours in the co- 
moving frame. For our present purposes, we determine the 
propagation of the pouch in a two step fashion. The first step 
determines the closed stream function contours in the resting 
frame of reference (as shown in Fig. 5, the so-called “stream- 
function gyre”) and calculates the corresponding zonal prop- 
agation of its area-centroid. The second step uses this first 
estimate of the propagation speed to determine the closed 
stream function in the the co-moving frame (also shown in 
Fig. 5) and obtain an updated propagation speed of the cor- 
responding centroid. In practice, the second step gives only 


a small correction (of a few percent) to the first guess for the 
propagation speed and for this reason we have neglected the 
second step for the calculation of the propagation speed in 
the majority of the experiments presented here. 

The advantage of the centroid method is that since the cen- 
troid is defined based on the spatial average within the closed 
streamfunction gyre in the resting frame (or pouch in the co- 
moving frame), high-frequency wind fluctuations exert little 
influence on the track of the centroid. As a result, the esti- 
mated zonal propagation speed (e.g.. Fig. 6) is slowly vary- 
ing with time and rather robust. Unless otherwise noted, the 
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□g. □ Stream function in the resting frame of reference (dotted 
contours) and in the co-moving frame of reference (solid contours). 
The thick solid line is the bounding contour of the wave pouch in 
the co-moving frame of reference, and the thick dashed line delin- 
eates the closed circulation in the resting frame (stream function 
gyre). The circle and the cross indicate the pouch center and the 
mass centroid of resting stream function gyre, respectively. 


(a) Longitude of the Centroid 



□ig. □ aQThe longitude of the mass centroid and LfoDthe propaga- 
tion speed (Cp) of the gyre-pouch based on its vorticity centroid 
as described in the text. The black curves represent the coarse- 
resolution simulation, and the green curves represent the warm-rain 
experiment. The abscissa is time (h). The ordinate is longitude in 
(a) and propagation speed (m s _ 1 ) in (b). 

time-varying propagation speed as obtained above is used to 
translate the frame of reference, and all pertinent diagnoses 
related to the vorticity dynamics and thermodynamics of the 
pouch region are carried out in the co-moving frame and are 
centered on the pouch centroid. 
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Comparison of pouch structure in the co-moving frame as 
initialized in the moist runs (Fig. 5) with that produced by the 
dry run (Fig. 2) indicate that both have a region of closed cy- 
clonic circulation centered at 18°N, with minor differences 
in kinematic structure near the periphery of the pouch. They 
share a common feature having a central “core” region of 
positive Okubo- Weiss parameter (following the sign conven- 
tion of DMW09) implying a region favorable for vortex ag- 
gregation and shape preservation of the meso-beta vorticity 
aggregate or “proto-storm.” 

□ numerical eC^eriments and supporting 
interpretations 

□ .1 DT □ revisited Dcoarse Resolution QEDkmD 

simulation 

In the control run, the relative humidity is specified as in 
KT81, so that the planetary boundary layer is uniformly 
moist with relative humidity values of approximately 85 , 
consistent with observations over the West Atlantic during 
the hurricane season. 

Figure shows time series for maximum wind speed 
at 850 hPa, surface pressure and 850 hPa relative vorticity. 
Three time series are shown for each of the three variables, 
corresponding to an increasing level of complexity in the rep- 
resentation of moist processes. The black curves denote the 
coarse resolution simulation using the Betts-Miller- anjic cu- 
mulus parameterization scheme; the green curves depict the 
high-resolution “warm-rain” simulation, and the blue curves 
depict the high-resolution simulation with full ice micro- 
physics (WRF Single-Moment 6-class scheme; Hong and 
Lim, 2006). The coarse resolution experiment will be dis- 
cussed here; discussion of the latter two simulations that ex- 
plicitly represent cloud processes will be deferred to the up- 
coming sections. 

The horizontal grid spacing in the coarse resolution sim- 
ulation is 28 km (or 0.25°) in contrast with the somewhat 
larger grid spacing of 0.625° in KT81. The maximum hor- 
izontal wind speed at 850 hPa attains minimal hurricane 
strength at approximately 1 44 h and decreases in intensity 
after 1 80 h. The minimum surface pressure drops slightly 
below 990 hPa soon after 1 44 h. The maximum 850 hPa rela- 
tive vorticity is 1 .5 x 1 0 -3 s _1 and is attained a little earlier at 
128 h. While the time series is broadly useful for summariz- 
ing the time evolution of certain fields, it does not reveal the 
cyclogenesis pathway. The central question to be addressed 
here is what is the transformational process between the ini- 
tial wave-like disturbance and the emergent tropical storm 

Figure 8 shows the evolution of stream function and rel- 
ative vorticity in the frame of reference that moves at the 
same speed as the easterly wave. At early times the ini- 
tially broad-scale cyclonic vortex within the wave’s sur- 
face low stimulates the formation of multiple mesoscale 
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□g. □ Time evolution of 850 hPa maximum wind speed (a), sur- 
face pressure (b), and 850 hPa maximum relative vorticity (unit: 
10 _4 s _l ) (c). The black curve is for the single-grid, coarse- 
resolution simulation, the green curve is the high-resolution, nested 
simulation with warm rain microphysics, and the blue curve is the 
high-resolution, nested simulation with ice microphysics. The ab- 
scissa is time (h). 


vortices that subsequently intensify in the lower troposphere 
by vortex tube stretching in regions of deep (parameter- 
ized) moist convection. Generally, these vortices outlive 
the convection that creates them and tend to merge with 
neighboring remnant vortices or become axisymmetrized 
(“scavenged”) by younger and more intense convective 
vortices. These simulated vortical structures generally have 
the strongest winds near the top of, or within, the (inflow) 
boundary layer. From a vorticity budget perspective (not 
shown here; see upcoming high resolution analysis for more) 
the mesoscale vortex enhancement process is such that the 
intensification of the near-surface horizontal winds leads that 
of the middle levels, and by this definition is always from the 
bottom up (Montgomery and Enagonio, 1998; Montgomery 
et al., 2006; Tory et al., 2006a, b and Tory et al., 2007). Fur- 
thermore, the mesoscale vorticity enhancement occurs most 
persistently near the intersection of the wave trough and the 
critical latitude for the moist waves (purple curve). This in- 
tersection point exhibits a striking “attractor- like” quality in 
the ensuing development process. 


A close up view of the evolving stream function (solid 
curves) and relative vertical vorticity (color) during the time 
interval 78 h— 1 1 1 h is highlighted in Fig. 8. The evolution 
from 84 h to 108 h shows an example of merging mesoscale 
vortices. In this experiment, deep convective activity and 
low-level vorticity enhancement continue more or less un- 
abated near the central latitude, resulting in a strong tropical 
storm strength vortex with maximum winds at 850 hPa of ap- 
proximately 35 ms' 1 after 144 h integration time (Fig. 7a). 
The moral of the story suggested by this experiment is that 
there is an upscale, “bottom up” development that occurs 
within the wave’s surface low, very near the critical latitude 
of the parent (dry) wave disturbance. In this example the 
critical latitude is located coincidentally where the vertical 
shear of the zonal flow is minimal (cf. Montgomery et al., 
2006) and the meridional gradient of the basic state vorticity 
changes sign. 

Figure 9a shows a snapshot of the Okubo- Weiss param- 
eter (defined in DMW09 as OW=f 2 — S 2 — S 2 ) at 96 h, a 
time at which a tropical storm-strength vortex has emerged 
from within the parent wave’s pouch. As shown by 
McWilliams (1984), OW>0 indicates that the flow is vor- 
ticity dominant, whereas OW<0 indicates that the flow 
is strain rate-dominant. Regions with large positive OW 
tend to be immune from the filamentary enstrophy cascade. 
In the initial stage of pouch evolution (not shown), the 
weak mesoscale vortices are characterized by weak, posi- 
tive OW embedded in a local environment which is strain 
rate-dominant (OW<0). As the mesoscale vortices become 
stronger, they exhibit a structure similar to a tropical storm: 
at this time there are three distinct VHTs forming a curved 
band near the center of the pouch and these are marked by 
regions of strong positive OW (red/yell ow/light-green). This 
positive OW region is surrounded by a strain-rate dominant 
region (purple/dark-blue). 

Figure 9b-d shows a snapshot of horizontal divergence 
and equivalent potential temperature (0 e ) at the 850 hPa level, 
and saturation fraction (column relative humidity between 
surface and 500 hPa) at 96 h. Whereas the 850 hPa diver- 
gence field is noisy, as is to be expected in a strongly convec- 
tive environment, strongly convergent regions (purple) pos- 
sessing magnitude on the order of 10 _4 s _1 are clearly ev- 
ident near the center of the pouch. These maxima are ap- 
proximately collocated with the intense cyclonic mesoscale 
vortices (cf. Fig. 8) and OW maxima at this time. The cor- 
responding field 0 e (Fig. 9c) has maximum values that are 
generally collocated with convergence maxima. Addition- 
ally, the 6 e field outside this region possesses a spiral band 
structure in which bands of elevated 0 e spiral cyclonically 
inward, lending support to the terminology “feeder bands” 
(Guinn and Schubert, 1993). Finally, the saturation frac- 
tion (SF, Fig. 9d) indicates a nearly saturated column of 
moist tropical air in the center of the pouch with maxima 
again collocated with strong convergence or updrafts. The 
“feeder bands” are more evident in the SF field, and they 
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□I1DQI] Evolution of stream function and relative vorticity (color) in the coarse-resolution (control) run from hour 78 to hour 1 1 1 in the 
co-moving frame of reference that moves at the zonal speed of the parent easterly wave. The unit of the vorticity is 10 -4 s -1 . The black 
curve is the wave trough axis and the purple line indicates the (linear) critical latitude (see text for details). The ordinate is latitude and the 
abscissa is longitude. 


are associated with parameterized cumulus convection (not 
shown). These features are consistent with T81 (their 
Figs. 5-8 and Fig. 22), who found multiple convective sys- 
tems within the surface low. 

In rotating moist convection problems that commence 
with a relatively dry middle troposphere, one expects to 
see negative 0 e anomalies emerge within and ust above the 
boundary layer in association with convective downdrafts 
that import (dry) low 9 e air into the boundary layer (e.g., 
Rotunno and Emanuel, 1987). However, with the etts- 
Miller- aniic cumulus parameterization scheme, downdrafts 


are not explicitly active so convective downdrafts are non- 
existent, and consequently, the G e deficits in the boundary 
layer are considerably weaker than their counterpart simu- 
lations with explicit downdrafts (e.g., M06). (The lack of 
explicit downdrafts in the T81 study has been a lingering 
concern by some tropical cyclone researchers. This restric- 
tion is removed in the upcoming high-resolution simulations 
in Sects. 3.2 and 3.3.) 

It proves useful to introduce a new diagnostic that high- 
lights the evolution of the maximum relative vorticity along 
each latitude near and within the wave’s cyclonic critical 
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Hour 96, 850 hPa 



332 333 334 335 336 337 336 339 340 70 75 80 85 90 95 98 

mill I I O W (a. unit: 10~ 8 s -2 ), divergence (b. unit: 10 -5 s -1 ), 0 e ( ) (c). saturation fraction from the surface to 500 mb (shading in □), and 
wind field (vectors in d). The ordinate is latitude and the abscissa is longitude. 


layer. Figure 10a shows cyclonic vorticity with values 
larger than 1.5x10 4 s -1 only. A concentrated and intense 
mesoscale vortex appears around 96 h (Day 4) near the origi- 
nal critical latitude of the parent wave (18° ). After this time 

the vorticity intensifies primarily in a highly localized spa- 
tial region near the critical latitude, and the emergent vortex 
starts a slow northwestward drift away from the original criti- 
cal latitude. The chaotic swarming of vorticity maxima about 
the critical latitude (prior to 96 h) followed by the highly fo- 
cused vorticity concentration is a striking feature. The multi- 
ple vortex activity exhibited prior to the attractor-type behav- 
ior is somewhat reminiscent of T. Fu ita’s infamous “suction 
vortices’ first hypothesized from retrospective Tornado dam- 
age surveys (Fu ita. Mystery of Severe Storms, 1991). The 
multiple vortices evident in this figure are obviously consid- 
erably larger in horizontal extent than tornado suction vor- 
tices, but these vortices are nonetheless smaller in horizontal 


scale than the tropical storm vortex and represent a class of 
analogous sub-vortex scale structures. This type of diagram 
will hereafter be referred to as a “Fu ita diagram” 1 . 

Figure 10b shows the Fu ita diagram of the analyzed satu- 
ration fraction (column relative humidity from the surface to 
500 mb). The emergent vortex, which develops into a trop- 
ical storm, forms in a region with high saturation fraction, 
efore storm formation the horizontal scale of the region of 

1 The link with Fu ita's work is in reference to the multiple vortex 
nature of the genesis process within the pouch rather than the mech- 
anism responsible for the multiple vortex structures. In the tornado, 
the origin of the suction vortices is thought to be due to a shear 
instability associated with a two-celled vortex. Prior to the genera- 
tion of a tropical storm, the small-scale cyclonic vorticity anomalies 
originate from convective instability and the associated vortex-tube 
stretching of like-sign vorticity of the local environment in the lower 
troposphere. 
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mm ii I I l “Fn ita Diagram:" (a) depicts the time evolution of the maximum relative vorticity along each latitude. The Jt-axis indicates the 
latitude of the vorticity maximum and the y-axis is time (from Day 0 to Day 10); (b) depicts the time evolution of the maximum saturation 
fraction from surface to 500 hPa. The ordinate is time (h). 
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n ull i I i i The 3-grid nested simulation: Relative vorticity and the corresponding (Earth-relative) streamlines at the initial time (a) and alter 
120 h integration (b). The two rectangles represent the two inner grids. The ordinate is latitude and the abscissa is longitude. 




high saturation fraction is comparable to that of the region of 
large cyclonic vorticity. After consolidation to a single mas- 
ter vortex there is a broader “funnel” of saturation fraction. 


□ □□□□[!] nu n ii i iirirn nn nmmtiiD (cmuxi) um 

□Dam GhODDD StTETlZBD 

As a simple way of testing both the robustness of the T81 
results and the new interpretation presented herein, we move 
progressively to a nested grid numerical simulation of the 
same experiment conducted above. The experimental set up 
is highlighted in Fig. 1 la. Three grid domains are indicated 
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by the solid black rectangles. The coarse horizontal grid 
spacing is 28 km (the same as before); the intermediate 
horizontal grid spacing is approximately 9.3 km; while the 
fine grid spacing is approximately 3.1 km. The horizon- 
tal extent of the outer grid is increased to two wavelengths 
(50° longitude), and the innermost grid contains only one 
full wavelength. For the time period chosen, the cyclonic 
trough region of the eastern wave spends the ma ority of its 
time within the intermediate resolution (9.3 km) grid domain 
where the etts-Miller- aniic cumulus parameterization is 
operative, and the cyclonic trough region of the western wave 
spends all of its time in the high resolution (3.1 km) grid 
with explicit moist convection operating therein (Fig. 1 lb). 
In this experiment, the WRF warm rain scheme is used, and 
the effects of ice microphysics are neglected. 

An advantage of this experiment is that one can assess 
the relative importance of explicit moist processes (warm 
rain and ensuing precipitation driven downdraft phenomenol- 
ogy) and parameterized moist processes within the same nu- 
merical experiment. For example, from the vorticity field 
shown in Fig. 1 lb, a concentrated tropical storm-like vor- 
tex has formed by 1 20 h on both the intermediate and fine 
grids. (The fine and intermediate grid fields have been in- 
terpolated up to the coarsest grid, 28 km.) Although the 
storms are located at approximately the same latitude, a con- 
solidated tropical-storm vortex looks more well developed in 
the intermediate grid. This is confirmed from the time se- 
ries of 850 hPa horizontal wind speed, minimum sea level 
pressure and 850 hPa relative vertical vorticity (not shown). 
The maximum horizontal wind speed and sea-level pressure 
deficit predicted on the fine grid are reduced relative to the 
coarse grid simulation. The primary reason for the delayed 
development in the warm rain scheme lies with the occur- 
rence of precipitation-induced downdrafts which import rel- 
atively low 0 e into the boundary layer, thereby interrupting 
the moistening and destabilization of the lower troposphere 
by sea-to-air fluxes of water vapor (e.g., Montgomery et al., 
2010). For the remaining discussion we present pertinent 
highlights of the warm rain simulation in as much as they 
support or refute the coarse resolution test and the proposed 
marsupial paradigm. 

□ □□□oamaDDDnn nrmrrnn mm nnn 
m i uni mu imaatHD 

As is evident in the translated stream function and precipita- 
tion rate patterns in Fig. 12, the initially broad-scale cyclonic 
vortex within the wave’s surface low stimulates the forma- 
tion of multiple precipitation events associated with the vor- 
tical hot plumes represented by the high-resolution grid (not 
shown). The horizontal purple line in each figure denotes the 
(linear) critical latitude defined by equating the zonal mean 
velocity to the zonal wave propagation speed (Fig. 6). At 
early times (^72 h) precipitation is widely distributed within 
the wave’s pouch. At later times, however, precipitation in- 


creases and becomes more and more concentrated near the 
center of the pouch. At long times (~ 1 20 h; not shown) max- 
imum precipitation is coincident with the pouch center and 
the intersection of the wave trough axis and the (linear) wave 
critical latitude. 

An example of the merger of the mesoscale vortices or 
their vortical remnants is shown in Fig. 13. At 82 h, there are 
two prominent cyclonic vortices in the northeast quadrant of 
the pouch. They move westward with the pouch. Mean- 
while, the weaker vortex (LI) rotates cyclonically around the 
stronger vortex (A), and it merges with the stronger vor- 
tex by 84.4 h. The resultant vortex has a stronger intensity 
and larger size than the individual vortices, and it is asso- 
ciated also with high surface 0 e as shown in the right col- 
umn in Fig. 13. This dominant vortex moves to the center 
of the pouch by 87 h (not shown). As confirmed by the 
Fu ita diagram shown in Fig. 14, the vortex development 
has transitioned from a swarm of approximately independent 
VHTs and their vortical remnants to a single master vorticity 
monolith immersed in a nearly saturated tropospheric col- 
umn (Fig. 14b). This attractor- like evolution within the wave 
pouch noted previously persists here. 

ased on the experiments presented in this paper, and 
in a companion study examining the cyclogenesis of pre- 
Hurricane Felix (2007) (Wang et al., 2010a, b), this attractor- 
like behavior appears to be a ubiquitous property of the simu- 
lated wave-to-storm transition process. The focal point is the 
nominal center of the elvin cat’s eye of the finite-amplitude 
tropical easterly wave (or the center of the wave pouch). 
From the perspective of linear wave dynamics, the prefer- 
ence for persistent convection near this focal point was al- 
ready anticipated theoretically, as well as observational ly us- 
ing the ASA/TRMM precipitation products, in DMW09. 
Recall that along the critical latitude the intrinsic frequency 
vanishes and fluid parcels are able to reside the longest there. 
During this time, parcels experience persistent moistening 
from the underlying ocean and latent heating associated with 
deep cumulus convection. ecause these fluid parcels are 
immersed in a cyclonic-rich vorticity environment, the latent 
heating will result in the stretching of the cyclonic vortex 
tubes passing though these parcels in the lower troposphere. 

arring unfavorable conditions, these two properties virtu- 
ally guarantee a rapid intensification of the local cyclonic 
vorticity near the center of the cat’s eye flow. The “sweet 
spot’’ of the cat’s eye (defined here as the intersection be- 
tween trough axis and critical latitude in the lower tropo- 
sphere) is thus the preferred location for tropical cyclogen- 
esis and the findings herein support this novel prediction of 
the marsupial paradigm. 

A mesoscale view of the wave-to-storm transition process 
is obtained next by calculating horizontal averages of key 
variables over a spatially limited portion of the wave pouch 


Atmos. Chem. Phys., 10, 10803-10827, 2010 


www.atmos-chem-phys.net/10/10803/2010/ 


181 


M. T. Montgomery et al.: Coarse, intermediate and high resolution numerical simulations 


10815 


Precip and 850 hPa Psi 






9 11 13 15 17 19 


OiED Evolution of translated stream function and precipitation (shading; unit: mm/h) from the inner grid (3 km resolution) simulation 
from Hour 72 to Hour 116. The zonal straight line indicates the critical latitude and the meridional black curve represents the trough axis. 
The ordinate is latitude and the abscissa is longitude. 


containing the sweet spot, and displaying the results in a 
time-height format. Figure 16 presents the results from such 
an analysis for a 2° latitude by 2° longitude box centered on 
the translating sweet spot. Shown are vertical distributions of 
area-averaged relative vorticity, horizontal divergence, rela- 
tive humidity and 0 e fields. Figure 15 highlights certain as- 
pects of this area-average analysis in a time series format. 

efore examining the details of Fig. 1 5, it is helpful to note 
that the fine grid is switched on at t 46 h and as a result there 
is some ad ustment of the fields on the fine grid until approxi- 
mately t 66 h. efore the fine grid is switched on, the area- 
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averaged cyclonic vorticity increases as in the coarse reso- 
lution simulation. Once the fine grid is switched on, how- 
ever, the deep convective processes trigger convective down- 
drafts (not shown), and for an interval of approximately 20 h, 
the area-averaged surface vorticity and G e weaken until the 
boundary layer has moistened sufficiently (cf. Fig. 4 of Ro- 
tunno and Emanuel, 1987). 

As shown in Fig. 1 5, after t 66 h, the near-surface vortic- 
ity (900 hPa eta) intensifies until t 90 h when the intensi- 
fication is interrupted briefly; amplification is then renewed 
with a somewhat larger rate. Also evident after t 66 h is 
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□llll I 1. 1 Evolution ofVHTs from Hour 82 to Hour 84.4 (3 km mesh). The left panels are 850 hPa relative vorticity (unit: 10 -4 s -1 ), the right 
panels are the surface 0 e ( ). and both are superimposed on the 850 hPa translated streamlines. The black curve is the wave trough axis, and 
the purple line represents the wave critical latitude. The two vortices discussed in the text are indicated by black arrows and letters (A and 
) in the left panels. The ordinate is latitude and the abscissa is longitude. 


that the amplification of near-surface vorticity leads the am- 
plification of mid-level vorticity (500 hPa eta). As in the 
coarse grid simulation, the simulated vortex development is 
thus from the bottom up. 

Once the fine grid is turned on the boundary layer RH 
always exceeds 85 (see Fig. 16c); during the grid transi- 
tion period the RH quickly increases to near saturation values 
(>95 ), remaining as such thereafter. Since the near-surface 
RH is nearly saturated after t 60 h, the near-surface RH is 
not useful for predicting the genesis time. It has been sug- 
gested recently that the increase in midlevel humidity may 
be a “trigger” for tropical cyclogenesis. The proposal is that 
the increase of midlevel RH generally precedes the amplifi- 
cation of surface vorticity and as result this might be a crucial 
observable for forecasting the timing of cyclogenesis ( olan. 


2007). It is clear from Fig. 15, however, that in the context 
of the tropical wave-to-cyclone transition problem as posed 
by T8 1 and by us, the increase of mid-level relative humid- 
ity (500 hPa RH) does not precede the amplification of near- 
surface vorticity (900 hPa eta); rather, it accompanies the 
spin up of near-surface vorticity, and both follow the com- 
mencement of the low-level convergence (900 hPa Divg). 
The evolution of relative humidity at 400 hPa or 600 hPa is 
similar to that at 500 hPa. 

The time-height evolution of divergence, relative vorticity, 
relative humidity and 0 e is shown in Fig. 1 6. The time-height 
plot of horizontal divergence shown in Fig. 16b depicts a 
fairly deep low- to mid-level column of divergence during 
the coarse-to-fine grid transition period (48 h< r<66 h); this 
is associated with the precipitation-induced downdrafts and 
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nrrm rm Fii ita diagram at high resolution: Time evolution of the 850 hPa maximum relative vorticity (left; unit: 10 -5 s _l ; derived from 
the 3-km horizontal grid output) and saturation fraction (SF) along each latitude. The x-axis indicates the latitude of the vorticity or SF 
maximum, and the y-axis is time (from Hour 25 to Hour 120). The thin curve in the left panel indicates the latitude where the zonal mean 
vertical shear between 200 hPa and 850 hPa is zero. 
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nmn rm Time evolution of 500 hPa relative humidity (unit: ), 

900 hPa divergence (unit: 1 0 -5 s~ 1 ), 500 hPa and 900 hPa relative 
vorticity (unit: 10 -4 s -1 ) averaged in a 2°x2° box following the 
moving pouch center from Hour 50 to Hour 120. 


the attendant low-level divergence patterns that accompany 
these downdrafts. Even during this transition time, how- 
ever, there is a clear surface convergence signature. Soon af- 
ter t 66 h, a moderately deep convergent column spanning 
the low- to mid-troposphere replaces the previous divergence 
pattern. The convergence within this layer intensifies slowly; 
the top of convergent layer also descends with time until the 
top is below the 850 hPa level for t >120h. Subsequently, 
all of the convergence is confined to a relatively thin bound- 
ary layer. (The upcoming section will suggest that this latter 
property is likely due to the absence of the ice microphysics 
in this simulation.) The divergence field at t 66 h is con- 
fined to above 450 hPa, but with time the level of zero diver- 
gence gradually descends to approximately 850 hPa at 120 h. 
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The time-height plot 9 e in Fig. 16 depicts a gradually de- 
creasing surface 0 e during the early phase of the genesis pro- 
cess (48 h <t <90 h) until a time at which rapid intensification 
of the emergent vorticity monolith is well underway. This is 
unlike the widely accepted view that the increase of the local- 
average vorticity should be accompanied by an increase in 
the local-average G e near the ocean surface, as expected in 
the familiar W1SHE mechanism of tropical cyclone intensi- 
fication articulated by Emanuel et al. (1994). This finding 
is supported by recent observational work (Molinari et al., 
2004) and theoretical/modeling work (Montgomery et al., 
2006; Tory et al., 2006b) indicating that tropical cyclogen- 
esis is not a WISHE process. These studies suggested that 
the amplification process during the genesis stage is domi- 
nated by deep and locally buoyant convective vortical plumes 
that are able to extract sufficient moisture from the boundary 
layer without significant augmentation of the sensible and 
latent heat fluxes associated with an increase in the surface 
wind speed on the vortex scale. 

□ □□□□ n ii li aim nj i iiti robDnmiD 

Similar to the procedure expounded in Haynes and McIn- 
tyre (1987) and the methodology outlined in Weisman and 
Davis (1998), Cram et al. (2002) and Davis and Trier (2002), 
we use the flux form of the vorticity equation in pressure co- 
ordinates and calculate area-averages of the vertical vorticity 
tendency equation in flux form: 

drj d dv d du 

— = — T-(ur)+co- G) — —(vt) — co— + F) (8) 

dt 3jc dp ay op 

where rj is the absolute vorticity defined as r]= f+dv/ dx - 
du/dy , and F and G represent the zonal and meridional 
components, respectively, of the effects of additional forces, 
such as friction and subgrid scale processes. Here p denotes 
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2°X2° Box Average (dOl) 




mill I 1 1 i The vertical distribution of the 2°x2° box averages of (a) relative vorticity (unit: 10" 5 s" 1 ), (b) divergence (unit: 10 -5 s _l ), (c) 
relative humidity (unit: ) and (p) 9 e (unit: ) following the moving pouch center (defined by the intersection of the local critical surface 

and the trough axis) from 1 lour 48 to Hour 120. The ordinate is pressure (hPa) and the abscissa is time (hour). Values below the surface are 
masked out between hours 87 to 120. 


pressure and it is understood that horizontal derivatives are 
taken on pressure surfaces. Since the sub-grid scale pro- 
cesses are typically more than an order of magnitude smaller 
than the friction terms associated with the surface boundary 
layer, we may estimate the frictional contribution to vortic- 
ity tendency by residual using the calculated tendency and 
advective terms. Consistent with the methodology used to 
create Fig. 16, the pertinent velocities used in Eq. (8) are de- 
fined relative to the nearly constant translation speed of the 
sweet spot defined in the previous section. To reduce the 
high-frequency fluctuations, a one-hour running time aver- 
age is applied to each budget term. 
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In addition to avoiding the inherent problem of measuring 
the small difference between large terms that can occur when 
using the material form of the vorticity equation (Haynes and 
McIntyre, 1987), the strength of this method lies in the fact 
that, through Gauss’s theorem, the tendency of the circula- 
tion for any enclosed area fixed in space can be written in 
terms of the line integral of the flux component normal to the 
boundary of the area (Davis and Trier, 2002): 


W = //^"=T v,_ "‘*V + "' f) ' sa ' <9) 

www.atmos-chem-phys.net/ 1 0/ 1 0803/20 1 0/ 


185 


M. T. Montgomery et al.: Coarse, intermediate and high resolution numerical simulations 


10819 


where C is the absolute circulation, A represents the area 
encompassed by the closed loop, 7 the horizontal velocity, 
and "/?=(— G, F). We will refer to the first term under the 
integral as the horizontal convergence of absolute vorticity 
(horizontal flux term). A comparison between the flux and 
the material forms of the vorticity equation (Holton, 1992) 
reveals that the horizontal flux term in Eq. (9) is related to 
the horizontal advection term plus the stretching term in the 
material form of the vorticity equation. The second term is 
an amalgamation of the vertical advection of vertical vortic- 
ity and the tilting of horizontal vorticity into vertical vortic- 
ity and will be hereafter denoted simply as the tilting-like 
(or non-advective) term (see Tory and Montgomery, 2008 for 
clarification of this terminology). For comparison, Davis and 
Trier (2002) refer to these terms as stretching and tilting, re- 
spectively. 

The budget terms are illustrated in Fig. 1 7. They are calcu- 
lated by averaging the vorticity flux divergences in a 2° x 2° 
box as done in Fig. 16, which by Stokes’ theorem is equiva- 
lent to the line integral of the corresponding flux through the 
lateral boundary. Figure 17a shows that the area-averaged 
vorticity tendency during the period between 66 h and 80 h 
is generally positive between the low- to mid-troposphere 
and the boundary layer, which is consistent with the weak 
convergence from the low- to mid-troposphere as shown in 
Fig. 16b. After t 80 h, however, Fig. 17b shows that the 
horizontal convergence of vertical vorticity becomes more 
focused in the lower troposphere and it is eventually largest 
in a shallow boundary layer (/ > lOOh) (the top of which 
we associate nominally with the depth of the strong conver- 
gence layer). Since the tilting-like term shown in Fig. 17c 
(which includes the resolved convective momentum trans- 
port) makes a second-order contribution to the calculated 
vorticity budget, and since the friction term (as calculated 
by residual in Fig. 17d) is maximized within this layer and is 
always negative there, the only process that can account for 
the spin up of cyclonic near-surface vorticity at these later 
times is horizontal convergence of vertical vorticity within 
the boundary layer. This is borne out clearly in Fig. 17b. 

The spin up of near-surface vorticity after t 1 00 h merits 
some discussion in light of recent studies arguing that fric- 
tion, when confined to the sub-cloud layer, should always 
cause a spin down of the near-surface flow (e.g., Raymond et 
al., 1998; Raymond and Sessions, 2007; Marni et al., 2008). 
While it is correct that friction generally retards fluid mo- 
tions, it will soon become clear that such reasoning is overly 
simplistic in vortical flows possessing VHTs whose aggre- 
gate effect is to establish a near-surface convergence of ab- 
solute angular momentum (Hendricks et al., 2004; Mont- 
gomery et al., 2006). For the present numerical simulation 
we have calculated the complete vorticity budget. In order to 
account for the spin up of the emergent tropical storm it was 
not necessary to assume a “boundary layer depth of 4 km” 
(cf. Raymond et al., 1998). 
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Recent work by Smith et al. (2009) and ui et al. (2009) 
have demonstrated that there are generally two spin up mech- 
anisms operating in tropical cyclones. The first is associ- 
ated with a diabatically-forced convergence of absolute an- 
gular momentum above the frictional boundary layer. The 
second is via a diabatically-forced but frictionally modified 
inertial inflow that is maximized within the boundary layer. 
In the latter case, the frictional boundary layer can gener- 
ate a positive spin up effect provided the convergence of ab- 
solute angular momentum acts fast enough to overcome the 
corresponding transfer of angular momentum into the under- 
lying ocean. This is essentially what happens post genesis 
(t > 100 h) in Figs. 1 6 and 1 7. 

The moral of the story is that the spin up of near-surface 
vorticity in the experiments conducted here results primarily 
from the near-surface convergence associated with the dia- 
batically forced inflow from the buoyant VHTs. Moreover, 
during the post-genesis/intensification phase, the leading- 
order surface spin up mechanism is from convergence within 
the frictional boundary layer, which is again driven primarily 
by the buoyant VHTs near the pouch center. This assessment 
is consistent with the findings reported in our companion pa- 
pers (Wang et al., 2010a, b) for the real-case simulation of 
the genesis of pre-Hurricane Felix (2007). 

nnnnnm nunm ii iiirnim mn nranrm(irmm) nm 
rannra rmnma n 

In this section we take up the issue of the modification of the 
above genesis sequence due to ice microphysical processes. 
The motivation for this section is manifold. One hypothesis, 
as extreme as it may seem, is that the results and interpreta- 
tion will change fundamentally when ice microphysical pro- 
cesses are included. To be sure, the use of rain-only physics 
provides an incomplete picture of the mesoscale processes 
occurring within the pouch; in particular, it does not address 
the hypothesized relative role of convective and stratiform 
convective processes in the genesis sequence (DM W09). We 
employ the same nested grid configuration used in Sect. 3.2 
and the same boundary layer and cumulus parameterization 
on the coarsest mesh (28 km). The only difference is that 

essler warm-rain microphysics scheme is replaced by the 
WRF single-moment, 6-class microphysics scheme (Hong et 
al., 2006) which includes ice, snow and graupel processes. 

The time series of the high resolution simulation with ice 
microphysics is summarized in Fig. 7. As in the warm rain 
experiment, the fine grid is activated at t 44 h. For all three 
metrics shown, the ensuing development is significantly de- 
layed relative to the warm rain experiment. It is not until 
approximately t 1 80 h that the ice microphysics simulation 
surpasses the intensity of the warm rain experiment, presum- 
ably because of the stronger downdrafts associated with ice 
processes (Srivastava, 1987). 
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□ HI Li 1 iVortieity budget (averaged over a 2-degree x 2-degree box following the pouch center) from Hour 48 to Hour 120. All fields are 
derived from the inner most grid (d03) simulation and a one-hour running mean was applied to reduce high-frequency fluctuations. The 
vorticity tendency unit is 10 -9 s -2 . The ordinate is pressure (hPa) and the abscissa is time (h). 


The mesoscale view of the wave-to-storm transition process 
is presented in Fig. 1 8, following the same methodology used 
to construct Fig. 16. Shown are vertical distributions of area- 
averaged relative vorticity, horizontal divergence, relative hu- 
midity and Q e fields. 

eginning around t 66 h, the near-surface vorticity 
slowly intensifies. This is accompanied by an amplification 
of the low- to middle troposphere vorticity field. As in the 
warm rain experiment, we see a deep convergent layer that 
becomes progressively shallower on the whole with time as 
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the system vortex intensifies. Consequently, the region of 
mean divergence deepens with time, but retains its maximum 
in the upper troposphere as expected. The convergence is 
maximized in the boundary layer and becomes pronounced 
when the system vortex begins to intensify rapidly (r > 1 44 h. 
Fig. 7). A time series of the near-surface and 500 hPa vortic- 
ity (not shown) shows that the amplification of surface vor- 
ticity leads the amplification of the 500 hPa vorticity. The 
development process is thus again from the bottom up. 
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2°X2° Box Average (dOl) 
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fi rm i i i S ame as Fig. 16, except for the simulation with ice microphysics. See text for details. 


Even with the standard WRF ice microphysics package, 
we do not see evidence to support the middle-level “humid- 
ity trigger” mechanism suggested by olan (2007). From 
Fig. 1 8c, we see again that in the context of the tropical wave- 
to-cyclone transition problem posed by T8 1 , the increase of 
mid-level relative humidity (500hPaRH) does not precede 
the amplification of near-surface vorticity (900 hPa zeta); 
rather, it accompanies the spin up of near-surface vorticity, 
and both follow the commencement of the low-level conver- 
gence (900 hPa Divg.) associated with the buoyant convec- 
tive towers. (Some further discussion of this issue is provided 
in our response to an Interactive Comment.) 


Finally, the time-height figure of 9 e (Fig. 18d) depicts an 
area-average surface 0 e that gradually decreases from t 66 h 
to approximately t 1 80 h; the latter is well into the construc- 
tion of the vorticity monolith. These results support the hy- 
pothesis that both the genesis and intensification process is 
not a WISHE process. The former is supported by Rotunno 
and Emanuel (1987) and the references cited earlier. The 
question of whether intensification can occur in an approxi- 
mately similar manner in the absence of the WISHE intensi- 
fication mechanism is answered in the affirmative elsewhere 
(Montgomery et al., 2010). 
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□□DU 1 1 iTime-radial evolution of the stratiform (left) and convec- 
tive (right) precipitation. The top panels show the area coverage by 
each rain type in percentage ( ), and the bottom panels show the 
azimuthal mean rain rate associated with each rain type. The ordi- 
nate is time (h) and the abscissa is radius with respect to the moving 
pouch center (km). 


□ □nnornrtriTn iTrTn nrm mria nirnTi 
nrTTTriTn mnnmrn 

We address next the matter concerning the relative role of the 
convective and stratiform convective processes in the genesis 
sequence. To do this we need a simple but physically-based 
method for identifying (deep) convection and stratiform pre- 
cipitation. Precipitation is categorized here as convective and 
stratiform following Tao et al. (1993) and raun et al. (2009), 
and the same algorithm is used in our companion study of 
Wang et al. (2010a, b). First, grid points with rain rate twice 
as large as the average of their nearest four neighbors are 
identified as convective cell cores. If a grid point is desig- 
nated as a convective cell core in this way, then its nearest 
neighbors (within one grid distance) are also designated as 
convective. Secondly, all grid points with surface rainfall 
rates greater than 25 mm h -1 are categorized as convective. 
To identify convective columns in which significant precip- 
itation is not reaching the surface, columns with maximum 
upward vertical motions larger than 5ms -1 or cloud liq- 
uid water below the melting level larger than 0.5 g kg -1 are 
designated also as convective. All remaining grid columns 


whose surface precipitation is greater than 0.1 mmlr 1 are 
categorized as stratiform; all grid columns with surface 
precipitation smaller than this rate are considered as non- 
precipitating. Our tests of this identification algorithm show 
that this method is sensitive to thresholds used to distin- 
guish the convective and stratiform rainfall. For example, if 
a sixteen-grid point average, instead of four-grid point aver- 
age, is used to identify convective cell cores, many more grid 
points will be categorized as convective. To avoid overesti- 
mating the convective contribution, some of the thresholds 
used in this study are higher than Tao et al. (1993) or raun 
et al. (2009). Our ob ective is to examine the quantitative 
evolution of stratiform and convective precipitation and their 
relative importance. We believe the algorithm used is ade- 
quate for this purpose. 

Figure 19 shows four plots. The two along the top depict 
a time-radius plot of the area percentage that is stratiform 
and convective within the wave’s pouch. The radius is de- 
fined relative to the center of the pouch and spans a domain 
of 300 km in extent. We see that the convective population 
is active throughout the pouch after t 48 h. The stratiform 
population on the other hand does not make a significant con- 
tribution until approximately t lOOh. In time, however, the 
convective population occupies a more substantial fraction 
of the area near the center of the pouch. This conclusion is 
confirmed in the second two figures along the bottom of the 
figure. The convective precipitation is maximized near the 
center of the pouch and progressively dominates the strati- 
form precipitation component. 

Figure 20 shows a snapshot of the rain rate and rain type 
distribution from the rain ice experiment at t 1 68 h. The 
top figure shows the rain rate in mm/h. The bottom figure 
shows the rain type. The heaviest rain rate is associated with 
the convective precipitation and it is spatially concentrated 
near the center of the pouch. 

These rain plus ice results point to a containment of per- 
sistent deep convection near the center of the pouch. This 
containment has a very important effect in the aggregate by 
producing mesoscale convergence at low-levels, which leads 
to low-level spin up. This behavior is unlike that of strati- 
form convection that selectively spins up the middle levels 
and spins down the low-levels (Tory and Montgomery, 2006; 
cf. Houze, 2004). 


□ □□□ □ a EDa □□ c EDc QD fill □ 

For tropical cyclogenesis in tropical easterly waves, the 
marsupial paradigm presented by DMW09 has begun to 
bridge the longstanding gap between observational ly resolv- 
able scales (synoptic, meso-a) and infrequently and insuffi- 
ciently observed meso-/3 and cloud scales (meso-y, VHTs) 
within and near the cyclonic trough of the parent easterly 
wave. sing a multiply-nested numerical model designed 
to test the key elements of the marsupial paradigm, we find 
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(b) Rain Type and 850 hPa Translated Streamlines (Hour 168) 
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Fig. 20. Snapshot of rain rate (a) and rain type distribution (b) from the ice microphysics run at Hour 168. In (b), green shading represents 
stratiform precipitation and red represents convective precipitation. 

broad agreement with the pioneering study of KT81 . More- 
over, the results suggest that the formation, aggregation and 
amplification of preferentially cyclonic vorticity anomalies 
by convective heating in the Kelvin cat’s eye of a critical 
layer constitutes the principal pathway for tropical cycloge- 
nesis in tropical easterly waves. Our findings directly sup- 
port the first hypothesis (H 1 ) as articulated in DM W09. This 
finding is demonstrated to transcend the particular horizontal 
resolution employed and is found to be robust to variations 
in cumulus parameterization, warm rain or ice microphysics 
formulations as represented in the standard “off-the-shelf” 
version of the WRF model. 

As discussed in DMW09, a key concept derived from the- 
ories of the nonlinear critical layer is that the cat’s eye ex- 
ists prior to organization and does not depend, structurally or 
otherwise, on the details of organization until TC intensifica- 
tion begins. On the other hand, as suggested in DMW09, 
the moist mesoscale processes may have feedback to the 
synoptic-scale wave and help to maintain and possibly en- 
hance the wave (H3). According to the second hypothe- 
sis (H2) of the marsupial paradigm, the formation of closed 
streamlines implies that the air within the materially closed 
space should be protected, to some degree, from lateral en- 
trainment of dry air. The closure serves not only to concen- 
trate vorticity by the aforementioned mechanism, but may 
account also for the accumulation of moist entropy and sat- 
uration of the tropospheric column by persistent deep moist 
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convection in the pouch - the same convection that facili- 
tates growth of cyclonic anomalies. This containment has a 
very important effect in allowing a convective type of conver- 
gence profile to dominate the closed region with a relatively 
smaller contribution from stratiform precipitation that is oth- 
erwise common in mesoscale convective systems (Tory and 
Montgomery, 2006; cf. Houze, 2004). 

We find that the maximum convergence in the central re- 
gion of the pouch occurs at a lower altitude than non-rotating 
organized tropical convection (e.g., Houze 2004), favoring 
rapid spin-up of one or more low-level vortices. The latter 
finding is consistent with Raymond and Sessions (2007) and 
moreover with the recent clarification of the spin up mecha- 
nisms operating in tropical cyclones by Smith et al. (2009). 
Using a more complete implementation of the WRF model 
that includes penetrative convection (hot towers, vortical or 
otherwise), ice microphysics and associated stratiform con- 
vection processes, the tendency for the deep convection to 
gradually dominate the stratiform convection near the wave’s 
sweet spot is verified here, consistent with H2. Our re- 
sults demonstrate that the marsupial sequence provides a di- 
rect “bottom-up” pathway to tropical cyclogenesis more or 
less independent of the mid-level transport of PV by typical 
mesoscale convective systems in the tropics (cf. Ritchie and 
Holland, 1997; Bister and Emanuel, 1997; Houze, 2004). 



(a) Rain Rate (mm/h) and 850 hPa Streamlines (Hour 168) 
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In our companion study (Wang et al., 2010a, b) we have 
conducted a numerical study of the genesis sequence of real- 
case Hurricane Felix that occurred in the main development 
region of the tropical Atlantic during the 2007 hurricane 
season. A full-physics implementation of the WRF model 
was employed in that study (including ice microphysical pro- 
cesses). Felix originated from an African easterly wave, and 
the low-level critical layer and its associated Kelvin cat’s 
eye formed from the interaction between the easterly wave 
and underlying Inter-Tropical-Convergence- one. As a re- 
gion of closed cyclonic circulation, the cat’s eye protected 
the moist air inside from dry air intrusion, and the formation 
of VHTs were favored in this vorticity- and moisture-rich en- 
vironment. Similar to the current study, the cat’s eye provides 
the focal point for the upscale aggregation of enclosed VHTs, 
and a tropical storm formed near the center of the cat’s eye 
via vorticity aggregation processes. 

In closing our two complimentary numerical studies, one 
idealized and one real-case example, provide strong support 
for the marsupial paradigm. Further observational, theoret- 
ical and numerical tests are ongoing and the results of this 
work will be reported in due course. In addition to these re- 
search lines, a new field experiment is forthcoming to test 
the theory in the Atlantic sector. The experiment and its re- 
lationship with the theory discussed herein is summarized in 
Appendix B. 


□ □[!□□□□ 
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The initialization method used here follows Kurihara and 
Bender (1980). The horizontal primitive equations in the o- 
coordinates may be written as (Phillips, 19 7) 


du o dp * dO 

dt ^ p dx dx 


(Al) 


dv 

d? +fu= - 


o dp* 
p dy 


00 

~dy 


(A2) 


where <7=-^, P* is the surface pressure, and 0 is the geopo- 
tential. Eliminating p using the equation of state, we have 


du 

~dt 


-fv = 


RT dp* 
P* dx 


00 

dx 


dv RT dp * 00 

— 4 - f u = — 

dt p * dy dy 


(A3) 

(A4) 


The geostrophic balance in the a-coordinates becomes 
fit x + /?TV//7/;* + V0 = O (A ) 


Assuming the zonal mean ow is in a state of geostrophic 
balance (for a at ocean surface) we then have 


dln/7* fu* 
dy ~ RT * 


(A6) 


00 

07 


-fu-RT 


d\np* 

dy 



(A7) 


where the subscript * denotes the value at the surface and the 
overbar denotes the zonal mean. iven p* and T(a) at a 
certain latitude 0o, and if T(o) is approximated by 7o(cr) in 
the above equations, we may solve for ~p * and 0. 

For an initial ow comprising a zonal ow plus a small, 
but finite amplitude easterly wave, the mass field must be ob- 
tained by solving the more general reverse balance equation 
in the a -coordinates 

V 2 <D + V(fl7V//ip,) = G (A8) 

where 


G = 2J(u,v) + fi;-up (A9) 

and where 7(w, The right-hand side of 

Eq. (A8) may be calculated from a given wind field. If, as 
above, T(o) is approximated by its vertical profile at the cen- 
ter latitude To(cr), then simplified Poisson-type equations for 
the surface pressure and geopotential are obtained for cr = l 
and o < 1 , respectively, as follows 

V 2 /np* = -T 1- (A10) 

R** o 

V 2 <D = G-7^g. (All) 

'*0 


For an initial horizontal wind field (given ~p * and 0 from 
the geostrophic inversion Eqs. A6, A7), the total fields (the 
zonal mean plus perturbation) may be obtained upon solv- 
ing Eqs. (A10, Al 1) by a relaxation method. Once 0 is de- 
termined, T is then obtained from the hydrostatic relation. 
The temperature field in Eqs. (A 10-All) can then be up- 
dated. This procedure is repeated several times until the tem- 
perature, surface pressure and geopotential fields no longer 
change appreciably. The fields so constructed thus solve the 
reverse balance equations. 


□□□crriDD 

□□□□□□□□□□ (20 □)) tiEmnmn 

Recent years have seen several field campaigns aimed at 
understanding the dynamics of tropical cyclone formation. 
These include the Tropical Cloud Systems and Processes 
(TCSP) in 200 (Halverson et al., 2007) and the ASA 
AMMA pro ect in 2006. Adding in the results of earlier ef- 
forts such as the Tropical Experiment in Mexico (TE ME ; 
Bister and Emanuel, 1997; Raymond et al., 1998) and even 
serendipitous observations of the early development of hurri- 
cane Ophelia in RAI E (Houze et al., 2006, 2009) and oc- 
casionally from reconnaissance aircraft (Reasor et al., 200 ), 
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and we have a collection of studies that have sampled pieces 
of a large and complex scientific puzzle. As discussed in 
DMW09, the puzzle begins with the formation of a tropical 
depression, a necessary meteorological precursor that results, 
in most cases, in the subsequent formation of a tropical cy- 
clone. 

So why should there be a new effort The first answer 
was given in the introduction to this paper the problem of 
tropical depression formation remains unanswered and is one 
of the great remaining mysteries of the tropical atmosphere. 
Perhaps the greatest shortcoming of previous campaigns is 
the limited sampling, both in space and in time. Therefore 
a second answer is that it is difficult to piece together snap- 
shots of tropical disturbances taken at different times. “ en- 
esis” (the formation of a tropical depression-strength vortex 
at sub-synoptic scales) often occurs between sampling times, 
or after disturbances move out of range. Third, there are very 
few observations of both precursors to genesis and the ensu- 
ing tropical cyclone formation process. The limited range of 
previous pro ects has meant also a limited phenomenological 
scope, regarding precursors to genesis and the multi-scale 
interactions needed for TC formation. We know now that 
synoptic-scale precursors take many forms (tropical waves, 
monsoon troughs, upper-tropospheric PV features, etc.). One 
of the sub-hypotheses for this experiment is that the genesis 
process itself should be relatively invariant on mesoscale and 
sub-synoptic scales in spite of these differences in synoptic 
precursor patterns. However, we do not know if this is cor- 
rect. Fourth, the marsupial paradigm of tropical cyclogene- 
sis discussed here embodies several new, unifying hypothe- 
ses that must be tested in order to confront the issues raised 
above. Finally, the combination of technological resources at 
our disposal is unprecedented. 

As discussed in DMW09 and in this paper, a key ingre- 
dient in the marsupial paradigm is the role of rotating deep 
moist convection and how convective organization in a rotat- 
ing environment differs from that of typical mesoscale con- 
vective systems in the tropics. Differences in cloud systems 
arising from the unique nature of the proto-vortex environ- 
ment imply a different pathway of dynamical evolution than 
what might otherwise occur in the absence of this environ- 
ment (e.g.. Figs. 16, 17, 18). The construction of a hurricane 
is a marvellous event, one that cannot be taken for granted 
in the current climate of Earth, nor in the perturbed climates 
to come. A unique cloud system structure and dynamical 
evolution, and a key role of precursor waves, contribute to 
this marvel in its early stages. The mesoscale aspects of this 
structure and evolution remain poorly understood. 

To address limitations of previous campaigns, the pro- 
posed PRE-Depression Investigation of Cloud-systems in the 
Tropics (PREDICT) aims to dramatically increase the spa- 
tial and temporal sampling of tropical disturbances prior to, 
and during, genesis. The PREDICT experiment is envisioned 
as a close collaboration between the Unites States ational 
Science Foundation ( SF), the ational Oceanic and Atmo- 
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spheric Administration ( OAA) and the aval Postgraduate 
School in Monterey, California. The primary research tools 
needed are the CAR -V aircraft, with altitude and range 
advantages over previous aircraft, and two OAA WP-3D 
aircraft each equipped with a tail Doppler radar, a lower fuse- 
lage (weather surveillance) radar, and the standard suite of 
ight-level and dropwindsonde observations that accompany 
normal hurricane and pre-hurricane reconnaissance. For the 
CAR -V aircraft, we propose to double-crew the aircraft 
for a portion of the field phase to allow sampling disturbances 
for as much as 16 out of 24 h. The pro ect will cover the 
ma ority of the Atlantic sector, including the Caribbean, and 
therefore will be poised to observe many forms of precur- 
sor disturbances and be positioned to uncover the common 
physical processes of genesis on the mesoscale. Finally, the 
pro ect aims to fully integrate the plethora of satellite data 
and derived products now available, coordinating the aircraft 
missions to maximize the total data coverage and physical 
interpretation. 
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